Theoretical phase diagram of ultrathin films of incipient ferroelectrics 
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We have used a first-principles-based scheme to compute the temperature-versus-misfit strain 
"Pertsev" phase diagram of ultrathin films of incipient ferroelectric KTaOa. Our results suggest 
that, at variance with the bulk material, KTaOs ultrathin films cannot be described as quantum 
paraelectrics. Rather, the behavior of the films is largely determined by surface/interface effects that 
favor ferroelectricity and the imperfect screening of the depolarizing fields. This leads to Pertsev 
phase diagrams that are qualitatively similar to those of normal ferroelectrics such as BaTiOa. 
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The need for miniaturized devices and the quest for 
knowledge in nanoscience have led to a flurry of activities 
in ferroelectric thin films that deepened our understand- 
ing of such low-dimcnsional systems (see Refs.[l|,[3,H and 
references therein). In particular, the crucial influence 
of surface/interface effects and mechanical and electri- 
cal boundary conditions on the properties of ferroelectric 
thin films is now well documented. For instance, recent 
studies revealed that the strain arising from the substrate 
can lead to ferroelectric phases absent in the bulk ^, ^, i6(| 
and that surface/interface effects can yield an asymmet- 
ric temperature-versus-misfit strain diagram [3, ^ . Simi- 
larly, residual depolarizing fields are now known to gener- 
ate anomalous effects in ferroelectric ultrathin films, such 
as the existence of stripe domains with remarkably small 
periods 

An important issue that remains to be fully understood 
concerns the properties of ferroelectric films made of the 
so-called incipient ferroelectrics, i.e., materials such as 
SrTiOs and KTaOa for which the bulk phase ferroelec- 
tricity is suppressed by quantum fluctuations 11, 12, T^ . 
The pioneering work in Ref. O revealed the existence of 
room-temperature ferroelectricity in strained 500 A-thick 
SrTiOa films, and can thus be taken to suggest that the 
coupling between strain and dipoles tends to prevail over 
quantum effects in nanostructures. One may thus won- 
der what the effect (if any) of quantum zero-point vibra- 
tions is in ultrathin films made of incipient ferroelectrics. 
In particular, the following questions are, to the best 
of our knowledge, currently unanswered: (1) Can quan- 
tum fiuctuations overcome the coupling between strain 
and dipoles in epitaxially strained films and, for exam- 
ple, suppress some phases? If so, which phases? (2) For 
zero misfit strain, are the films paraelectric or ferroelec- 
tric? (3) What does the coexistence of quantum and 
surface/interface effects lead to? (4) What are the con- 
sequences and signatures of the zero-point motion of the 
ions in a fllm experiencing a depolarizing fleld? 



The aim of this Letter is to answer these questions by 
determining and analyzing the temperature-versus-misfit 
strain phase diagrams of KTaOs ultrathin films under 
different electrical boundary conditions, using effective 
models of the films derived from first-principles. As we 
will see, many interesting features will be revealed and 
explained. 

First-principles-based effective Hamiltonian methods 
have been very fruitfully employed to study ferroelec- 
tric thin films before (see, e.g., Refs. i, i, [ia, andlH). 
Here we used such an approach to study KTaOa ultrathin 
films grown along the [001] pseudo-cubic direction (cho- 
sen to be along the z-axis), having K-0 terminated sur- 
faces/interfaces, and being 28 A thick. Technically, the 
films are modeled by 10x10x7 supercells that are peri- 
odic along the x- and y-axes (which lie along the [100] and 
[010] pseudo-cubic directions, respectively) and contain 
seven Ta02 (001) layers stacked along the non-periodic 
z-axis. The total energy of such a supercell is written as 



£tot = &cff (ui, Vi, 1]) + /3Ed • ^ Z*Ui , 



(1) 



where Ui is the local soft mode in the unit cell i of the 
film, such that Z*Ui yields the local electrical dipole, and 
/3Ed is a term related to screening of the depolarization 
field, as discussed below. The {vi}'s variables describe 
inhomogeneous strains within the supercell [13] . ij is the 
homogeneous strain tensor, which is particularly relevant 
to mechanical boundary conditions since epitaxial (001) 
films are associated with the freezing of some components 
of 7] (in Voigt notation), i.e., ?76 = and ?7i = 772 = S, with 
d being the value forcing the film to adopt the in-plane 
lattice constant of the substrate [1, @, 0]. In practice, 
(5=(asub — aKTa)/aKTa, whcrc Osub is the in-plane lattice 
parameter of the substrate and axTa — 3.983 A is the K 
cubic lattice constant of bulk KTaOa used in Ref. [12. 

The expression for fnoff, the intrinsic effective- 
Hamiltonian energy of the film, is given in Ref. in. The 
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first-principles-derived parameters for bulk KTaOa are 



those of Ref. [1^, except that here we use the formula 



for the dipole-dipole interactions for thin films under 
ideal open-circuit (OC) boundary conditions derived in 
Refs.|l3and[ll Such electrical boundary conditions nat- 
urally lead to the existence of a maximum depolarizing 
field inside the film (denoted by Ed) when the dipoles 
all point along the [001] direction. The second term of 
Eq. (1) mimics a screening of Ed via the parameter [3. 
More precisely, the residual depolarizing field resulting 
from the combination of the first and second terms of 
Eq. (1) has a magnitude equal to (1 — /3) Ed- As a result, 
/3 = corresponds to ideal OC conditions, an increase (3 
lowers the magnitude of the resulting depolarizing field, 
and /3 = 1 corresponds to ideal short-circuit (SC) condi- 
tions for which the depolarizing field vanishes. Note that 
Ed depends on the dipole configuration, is exactly de- 
rived at an atomistic level (following the procedure intro- 
duced in Ref. llTl). and is self-consistently updated during 
the simulations. 

We should note here the main approximations made 
in the construction of the above described model for 
ultrathin films of KTaOa under varying electrical and 
mechanical boundary conditions. Firstly, as in previous 
works [1, [3] , our f tot relies on a simple truncation, at 
the surface/interface layers, of the interactions existing 
in bulk KTaOa. While this approximation is admittedly 
crude, the resulting ftot mimics in a qualitatively cor- 
rect way the enhancement of the surface/interface po- 
lar modes that occurs in the most promising compounds 
(see, e.g^ the experimental work on Bao.sSro.sTiOs films 
in Ref. |7j or the ah initio results for PbTiOs films sand- 
wiched by Pt electrodes in Ref-IH). In our simple model, 
such an enhancement happens because our truncation re- 
moves short-range interactions that oppose the onset of 
such local dipoles. At a quantitative level, this approx- 
imation may well be underestimating the enhancement 
for some choices of electrodes and overestimating it for 
others. At any rate, as the interesting effects resulting 
from our simulations are very pronounced, we are confi- 
dent they can be taken as reliable qualitative predictions. 

Secondly, we made calculations for ideal electrodes 
leading to SC boundary conditons (i.e., /3= 1), and also 
for realistic situations for which an appropriate different 
value of (3 must be chosen. Let us mention two pieces 
of information to justify such a choice: (i) the ab intio 
studies of Sai efc al. [l^ , who found a 97 % cancellation of 
the depolarizing field in some systems, and (ii) ongoing 
unpublished work of one of our co-authors (L. BcUaiche) 
who, by comparing measurement and simulations, has 
estimated that some experimental set-ups correspond to 
a screening of about 98 % of the maximum depolariz- 
ing field. Here, we made a conservative choice and used 
f3 = 0.96. We thus expect the striking effects predicted 
will be robust as far as this approximation is concerned. 

As in Ref. [13, Stot is used in two different kinds of 



Monte-Carlo (MC) simulations: classical Monte Carlo 
(CMC) jl^, which does not take into account zero- 
point vibrations, and path-integral quantum Monte Carlo 



(PI-QMC) [21|, 122 



I, which includes the quantum- 
mechanical zero-point motions. Consequently, compar- 



ing the results of these two different Monte-Carlo tech- 
niques allows a precise determination of quantum effects 
on properties of the considered KTaOa ultrathin films. 
We typically used 30,000 MC sweeps to thermalize the 
system and 70,000 more to compute averages, except at 
low temperature in PI-QMC where more statistics are 
needed. Further details about the PI-QMC technique 
may be found in Refs.IHIil, and l23l. 
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FIG. 1: Temperature versus misfit strain diagram of a 28 A- 
thick KTaOs ultrathin film under ideal SC conditions (/3 = 1). 
Solid and open symbols refer to CMC and PI-QMC simula- 
tions, respectively. Serving as guides for the eyes, the solid 
and dashed lines are linear fits of the CMC and PI-QMC data, 
respectively. The inset shows a magnified view of the region 
in which the r phase occurs. 

Figure 1 displays the "Pertsev" (that is, the tempera- 
ture versus misfit strain S) phase diagram for the 28 A- 
thick KTaOs thin film under ideal SC conditions, as re- 
sulting from CMC (solid symbols) and PI-QMC (open 
symbols) simulations. One can clearly see that both 
kinds of simulations generate the four phases that also 
appear in the Pertsev diagram of (001) BaTiOa ultrathin 
films 0, 0] . These phases are: a paraelectric p state at 
high temperatures; a ferroelectric tetragonal c phase for 
intermediate temperatures and compressive and weakly 
tensile strains, having a polarization lying along the [001] 
growth direction; a ferroelectric orthorhombic aa phase 
for low temperatures and large tensile strains, in which 
the polarization is parallel to the in-plane [110] direction; 
and a ferroelectric monoclinic r phase for the lowest tem- 
peratures and intermediate tensile strains, for which the 
polarization direction continuously rotates from [001] to 
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[110] as S increases. 

Figure 1 further displays four remarkable features. (1) 
Quantum effects begin to appear at temperatures below 
wlOOK (we numerically found that the transition lines 
predicted by CMC and PI-QMC overlap above 100 K), as 
in buUk KTaOa [l2|- As a result, and as emphasized by 
the inset of the figure, the c-to-r and, to a lesser extent, 
the aa-to-r phase boundaries are the only boundaries 
affected by quantum fluctuations (PI-QMC simulations 
performed at temperatures higher than 100 K, at 0.5% 
and 1.3% misfit tensile strains, are indicated with arrows 
in Fig. 1, and reveal that the high-temperature p-to-c 
and p-to-aa phase transitions are unaffected by quantum 
fluctuations). The effects are, in any case, minor. Essen- 
tially, there is a small reduction of the strain range in 
which the r phase occurs: in the low-temperature limit, 
it passes from about -1-0.4 %< S < -1-1.3% at a classi- 
cal level to about -f0.5%< S < -1-1.3% when quantum 
fluctuations are included. (2) As in BaTiOa ultrathin 
films 0, the p, c, aa, and r phases meet at a sin- 
gle four-phase point, which occurs here at a tempera- 
ture of about 160 K and a tensile strain of about 1.12% 
for both CMC and PI-QMC simulations 0. (3) The 



phase diagram is asymmetric with respect to zero misfit 
strain. Such an asymmetry is hinted in experiments on 
Bao.sSro.sTiOa films and is related to the enhance- 
ment of the z-component of the local dipoles at the sur- 
faces/interfaces. As mentioned above, in our simulations 
the enhancement occurs because the polar local modes at 
the surfaces/interfaces are partl y fr ee from energetically- 



costly short-range interactions [25|. (4) For zero misfit 
strain, the paraelectric-to-c transition occurs at about 
460 K for both CMC and PI-QMC simulations. 

The high transition temperature (Tc) obtained for the 
film at zero misfit strain starkly contrasts with the be- 
havior the bulk material, for which the paraelectric-to- 
ferroelectric transition occurs at about 30 K at a clas- 
sical level and vanishes when quantum effects are con- 
sidered [13] ■ The physical origin of such a high Tq lies 
on the above mentioned dipole enhancement at the sur- 
face/interface of the films. The striking consequence of 
such an enhancement is that these ultrathin films are 
no longer incipient ferroelectrics, but display a Pertsev 
phase diagram that is in essence identical to that of a 
normal ferroelectric like BaTiOs. 

It is also interesting to realize that SrTiOa has a lat- 
tice constant about 2% smaller than KTaOa, and that 
a linear extrapolation of the data in Fig. 1 suggests that 
growing KTaOa ultrathin films on a SrTiOa substrate 
(for which S = —2 %) should lead to a Tq that could be 
as high as 1000 K. In other words, our simulations sug- 
gest that growing an ultrathin film made of an incipient 
material (that is, KTaOs) on a substrate made from an- 
other incipient ferroelectric (namely, SrTiOa) may result 
in a ferroelectric compound with a high Tq and, thus, 
a large polarization along the growth direction at low 



temperatures. This would be the result of, not only sur- 
face/interface effects, but also the well-known coupling 
between dipoles and compressive strain 
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FIG. 2: Temperature versus misfit strain diagram of a 28 A- 
thick ultrathin film KTaOa under a residual depolarizing field 
(/3=0.96). Solid and open symbols refer to CMC and PI- 
QMC simulations, respectively. The solid and dashed lines 
are linear fits of the CMC and PI-QMC data, respectively, 
and are guides for the eyes. The inset is a magnification 
of the region in which the r phase occurs, with the PI-QMC 
data point for the c-to-r transition temperature at zero misfit 
strain being a guess resulting from the interpolation of the PI- 
QMC simulations for the lowest computed temperatures. 

We now discuss how depolarizing fields, always present 
in reality as, e.g., the metallic electrodes are never ideal 
[2^, affect the Pertsev phase diagram of KTaOa ultra- 
thin films. As explained above, we have chosen f3 = 0.96 
in Eq. (1), i.e., a situation that corresponds to a screening 
of 96 % of the maximum depolarizing field. The resulting 
Pertsev diagram is shown in Fig. 2, for both CMC and 
PI-QMC simulations. By comparing the CMC results 
from Figs. 1 and 2, it is clear that the depolarizing field 
tends to suppress the z-component of the polarization, as 
expected. Indeed, this effect dramatically decreases the 
p-to-c transition temperatures; for instance, the transi- 
tion occurs around 60 K at zero misfit strain, i.e., the Tq 
is 400 K lower than the one obtained for ideal SC. It also 
considerably extends the region of small tensile strain 
associated with the aa phase, which exhibits an in-plane 
polarization that does not directly feel the depolarizing 
field, and significantly narrows the r-phase region. 

As a result, the phase diagram becomes more sym- 
metric around the zero misfit strain and the four-phase 
point of Fig. 1 shifts towards lower temperature (about 
50 K for CMC simulations [li]). As shown in Fig. 2, such 
a small transition temperature implies that the quantum 
fiuctuations (which are appreciable only at temperatures 
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below 120 K, and become more pronounced as the tem- 
perature decreases) have two remarkable effects that do 
not occur in the film under SC conditions. First of all, 
the zero-point vibrations now affect the p-to-c and p- 
to-aa phase boundaries (for small strain) - in addition 
to the boundaries involving the r phase. Secondly, quan- 
tum fluctuations significantly reduce the stability range 
of the r-phase. This suggests that, for slightly worse elec- 
trodes (i.e., for (3 < 0.96), quantum effects could well lead 
to the disappearance of the four-phase point in favor of a 
f/iree-phase point in which "only" the p, c, and aa meet. 

In summary, we have used a first-principles-based 
scheme to determine the temperature- versus-misfit strain 
"Pertsev" phase diagram of ultrathin films of incipi- 
ent ferroelectric KTaOa under short-circuit-like bound- 
ary conditions. We have performed both classical and 
quantum Monte Carlo simulations, the ionic quantum 
fluctuations being fully taken into account in the latter. 
This has allowed us to investigate the competition of the 
various effects present in the films, namely, the misfit 
strains and surface/interface effects that arc known to 
favor the occurrence of ferroelectric phases, and the zero- 
point vibrations and residual depolarizing fields that tend 
to suppress ferroelectricity. 

We have found that, for ideal electrodes, quantum fluc- 
tuations have a negligible effect; the phase diagram is 
essentially determined by the other factors mentioned 
above and closely resembles what is found in, e.g., 
BaTiOa films. Our results thus suggest that, in ultrathin 
film form, quantum paraelectrics such as KTaOs behave 
qualitatively in the same way as strong ferroelectrics such 
as BaTiOa. We have also simulated realistic (thus imper- 
fect) electrodes and found that the corresponding depo- 
larizing fields have a great impact on the calculated phase 
diagram; in particular, the transition temperatures of the 
c and r phases are significantly reduced. As a result, 
quantum mechanical effects become more important in 
this case, and may alter the phase diagram qualitatively, 
e.g., by suppressing the r phase. 

We believe the present study leads to a better under- 
standing of ferroelectric thin films, and hope it will stim- 
ulate experimental work aimed at confirming our predic- 
tions. 
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